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Abstract 

We investigate a production process e + p — ► tX — > bW + X at HERA, where 
we consider a decay mode b — > v e d of the sbottom in the framework of an 
R-parity breaking supersymmetric standard model. Both processes of the stop 
production e + d — > t and the sbottom decay b — > v e d are originated from 
an R-parity breaking superpotential X' l ^ l L\Q^D c i. One of signatures of the 
process should be a large missing transverse momentum plus multijet events 
corresponding to hadronic decays of the W . It is shown that the signal could 
appear as an event excess in the charged current (CC) processes e + p — > vX 
with the high Q 2 at HERA. We compare expected event distributions with the 
CC data recently reported by the HI and ZEUS groups at HERA. Methods 
for extracting the signal from the standard CC processes are also discussed. 
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1 Introduction 



Recently reported |l], |2| event excess at the high Q 2 region in the deep inelastic (DIS) 
process e + p — ► e + qX has very much excited theorists. Among its theoretical consid- 
erations proposed so far the supersymmetric (SUSY) model based approach seems to 
be an appealing possibility. It could be interpreted by the scalar top (stop) production 
e+p — ► tX — » e+dX Q in SUSY models with the R-parity breaking (RB) interactions 
proportional to A' 131 [Q, [|, [6|, [?], ^, |9| []. The excess events seem to be broadly distributed 
around 200GeV ~ 250GeV in the mass (M = ^/xs, where x is the Bjorken parameter). 
In a previous paper M we have shown that it could be simulated by our specific scenario 
if we consider almost degenerate two mass eigenstates t\ and ti of the stops with masses 
m~ ~ 200GeV and m~ ~ 230GeV. 

tl 52 



DO |TT[ and CDF [12] groups have recently obtained preliminary bounds on a lepto- 
quark mass as M LQ > 194GeV and 210GeV for Br(LQ e ± g) = 1.0, respectively. These 
new results suggest that the stop (or leptoquark) with Br(t — > e + d) = 1.0 and m~ ~ 
200GeVJias been excluded. Consequently stops could decay into not only e+d but also 
bW\ or bW via R-parity conserving interactions as far as there exists a chargino W\ or a 
scalar bottom (sbottom) b lighter than the stops. So we can reconcile the interpretation 
of the high Q 2 anomaly at HERA data in terms of stop with the Tevatron bound provided 
that Bi(t^ e+d) = 1- Br(t -»• bWi or bW) <, 0.8. 

Taking into account such situation we are naturally led to examine the processes (1) 
e+p — * bW\X or (2) e+p — > bWX. They will presumably give us unique experimental 
signatures of our scenario. We have shown that a possible typical signature of (1) would 
be a high Pt charged lepton plus jet(s) with a large missing transverse momentum fx, 
i.e., one of the signals to be detected at HERA is characterized by the high Pt spectrum 



of muons [14, 15]. Note that the lightest neutralino Z\ and b possibly decay into i?-even 
particles via only non-zero RB coupling A' 131 . Altarelli et al. have pointed out that another 
possible signature of (1) could be a large IpT plus mutijets, where the chargino decays into 
a neutrino plus quarks |jl3|| . Moreover, they have shown that such a sigal could be observed 
as an event excess at the high Q 2 region in the charged current (CC) DIS, e + p — > vqX : 
for a specific sparticle mass spectrum m~ ^ 70GeV and m~ ~ 180GeV. 

In the present letter we investigate the process (2) mentioned above, i.e., we are con- 
cerned with a case where the sbottom b is lighter than the stops and has a mass smaller 
than about 120GeV. Here we assume sufficiently heavy chargino W\ (<; 200GeV) and neu- 
tralino Z\ (^ lOOGeV) for simplicity. In this case the stops cannot decay into bW\ and the 
sbottom can only decay into v e d via the RB coupling A' 131 . First we will show that such a 
light sbottom has not been excluded from the present experiments. One of characteristic 
signatures of the process (2) would be a large missing transverse momentum plus multijets 
corresponding to hadronic decays of the W. It could be observed as an event excess at 
the high Q 2 region in the CC DIS. We will compare expected event distributions with the 



CC data recently obtained by the HI ||, lli and ZEUS yj] groups at HERA. 



* Another interesting possibily for "strange stop" scenario, e + + s — > with non-zero coupling X' Vc 
discussed in detail by J. Ellis et al. |10|. 
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2 Model 



We are based on the minimal SUSY standard model (MSSM) with an RB superpotential 

W H = \' 131 L 1 Q 3 Dc l , (1) 

where 1 and 3 are generation indices. We can immediately obtain Yukawa-type interactions 
from the superpotential ([[]) , 

L = A' 131 (tL~dP L e + r L dP L t + dR^P L t - V L dP L v e - v e dP L b - <T R l4P L b) + h.c, (2) 

where Pl reads the left handed chiral projection operator. The first term in the Lagrangian 
@ will be most suitable for the ep collider experiments at HERA because the stop tL will 
be produced in the s-channel for e + -d sub-processes. The fourth term will contribute to 
the sbottom decay, &l — > v e d. Note that the stop it, cannot couple to any neutrinos via 
the RB interactions. That is, no excess is expected in the standard CC events ep — * uqX 
in the squark scenarios. As has been pointed out, however, the CC like events could be 
expected if we consider some appropriate decay chains to the stops (see [O] and later). 



Here we pay attention to a fact that the stops [sbottoms] (tL, £r) [(&l, are naturally 
mixed each other fl9j| due to a large Yukawa coupling of their partner quark and the mass 
eigenstates (t±, t 2 ) [{bi, b 2 )] are parametrized by a mixing angle 9 t [Ob], 

tL = hcosOt -t 2 sm6 t , (3) 
&l = &i cos 9b — 62 sin Ob . (4) 

The interaction Lagrangian (||) can easily be rewritten in terms of the mass eigenstates. In 
particular, we should note that both stops t\ and t 2 could be produced by the eq scattering 
through 

Lj ed = A' 131 (cos OthdPLe - sin9 t t2dP L e) + h.c. (5) 
The most stringent upper bound on the coupling constant A' 131 comes from the atomic 



parity violation (APV) experiments [18, [2(j, 21, 22, 23, g4|. In the following discussion we 



adopt rather modest limit, A' 131 ^ 0.1. We will come back to this point later. 



3 Light sbottom scenario 

Interestingly, in the MSSM there exists a theoretical upper bound on the mass of the 
lighter sbottom 61. Since the left handed stop tL and sbottom 6l form a SU(2) doublet, 
their masses include the same contribution from a soft scalar breaking mass mQ 3 , 

= ffiQ 3 + mf + m 2 z cos 2(3 ^ — e u sin 2 Qy/^j (6) 
= cos 2 9 t m~ + sin 2 9tm~ , (7) 



TO?- = TO, 



q 3 + ml + m 2 z cos 2 p(^-^-e d sin 2 9 (8) 



cos 9bm~ + sin 9bm~ , (9) 
61 62 
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where the last terms in (6) and (8) stand for the D-term contributions. Combining for- 
mulae (6) through (9), we obtain an upper bound on the sbottom mass to 2 - as, 

b\ 



< 



< 



cos 2 9tm~ 



ml 

12 



nit 



sin 2 9trrv~ 

*2 



to 2 



m t 



mt\ cos 23. 



ml 



miY cos 23 



(10) 
(11) 



Thus the lighter sbottom b\ cannot be heavy for relatively light t 2 . For example, b\ will be 
lighter than t\ (m~ ~ 200GeV) for to^ ^ 250GeV even in the extreme case of tan 3 = 00 
(cos2/? = -l). 

Next we examine the decay modes of the stop. In the MSSM, the stop could be 
lighter than the other squarks in the first and second generations and the gluino. It can 
decay into various final states : t — ► tZ^ (a), bWi (b), Wb (c), blv (d), bvl (e), bW 
(f), bfJZ k (g), cZi (h) and ed (i), where Z k {k = 1 ~ 4), Wi{i = 1,2), v and J, 
respectively, denote the neutralino, the chargino, the sneutrino and the charged slepton. 
(a) ~ (h) are the .R-parity conserving decay modes, while (i) is only realized through the 
RB couplings (||). If we consider the RB coupling with A' 131 > 0.01, the decay modes (d) 
to (h) are negligible due to their large power of a arising from multiparticle final state or 
one loop contribution. Moreover, in the present case (mj ~ 200GeV) the mode (a) will 
kinematically be suppressed. Then only two body decay modes (b), (c) and (i) are left for 
our purpose. 

Here we assume the stops can decay through the modes (c) and (i). In other words, we 



consider sufficiently heavy charginos, 



/?? 



Wi 



> to 



to& i> 220GeV, and a light sbottom, 



mr < mr 
01 ti 



m\v ^ 120GeV. The former condition can be easily realized if we take such 
large soft breaking gaugino mass as M2 ^ 250GeV (and \fi\ ^> lOOGeV). As a result of 
this assumption, the lightest neutralino Z\ becomes inevitably heavier than about HOGeV 
because m 



Z-i. 



0.5m^ in the MSSM. Consequently, an R-parity conserving decay mode 



of the light sbottom, b\ — > bZ\, is almost kinematically forbidden. We can naturally 
consider b\ — > D e d as the dominant decay mode. It should be noted that the assumption 
of relatively light sbottom is consistent with the theoretical upper bound on mass (|ll|) for 
the light t 2 like m~ ~ 230GeV. 

Decay widths of U — > bjW + can be expressed as 



9 , f I 2 Pj 
levr 17 ^ 1 ml 



mw 



2 2 
m — h m~ 



+ 



2 

mA 

U 



mA 



tn 



w 



(12) 



where pj 



2 

mA 



E: 



(mA + mA 
v k b. 



m 



and (fn, f 12 , /21, f22) 



(+ cos 9t cos 9b, — cos 0t sin — sin t cos +sin0 t sm9b). Note that T(ti — > b\W + 
is proportional to cos 2 9^. In Fig.l, we show the branching ratio Br(ti — > e + d) = 1- 
Br(t — > b\W + ) as a function of m~ , where we take A' 131 = 0.1, m? 



210GeV and 9 t 



1.0. We can find that the desirable branching ratio Br(ti — > e + d) ^ 0.8 can be obtained 
for a wide range of the mass m~^ ^ 120GeV. It depends sensitively on the mixing angle 
9b through a factor cos 2 9b- 
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It is natural for us to ask whether or not such a light sbottom, m~ 5s 120GeV, has 
already excluded experimentally. In fact, the CDF and DO groups at Tevatron have 
excluded squarks lighter than about 200GeV [22]. We should note, however, that the 
limit has been obtained on the basis of assumptions of five degenerate squarks with tu~ l 
= fri~ R . That is, the expected cross sections will be significantly reduced by a factor ten 

when one sbottom b\ is lighter than the other squarks. At present, the most stringent 
mass bound on the sbottom with the dominant decay mode b\ — > v e d is 

m~ ^ 90GeV. (13) 

This has been obtained from a stop mass bound m~ ^ 90GeV for m~ = through search 



for a process, pp — > tit±X — > ccZ\Z\X at the DO 25]. Its physics background will be (1) 
the production processes for the stop and sbottom are almost the same at the Tevatron, 
to which SUSY-QCD diagrams are mainly contributed, and (2) experimental signatures 
of the stop and sbottom are the same, i.e., large fx plus jets. 

Next we should examine constraints from the precision measurements at LEP1. Po- 
tentially, contributions from the ti and bi to the Ap could become large if their masses 



are not so different from the weak mass scale of mz [||, 1C]. We have checked, however, 
contributions from the ti and bj to the Ap could become small to such extent as 1 x 10 -3 
even for m~ = lOOGeV, mj^ = 200GeV and mj 2 = 230GeV if we take large mixing angles 
Bt — 6b <^ 1-0 and a large tan/3 ^ 10. Here we have used the formula for Ap including 
mixings of both stops and the sbottoms pf| . 



4 Total cross sections 



The analytical expression for the differential cross section of the sub-process e + d — » b\W + 
is written as follows, 



d± _ g 2 A' 13 i 2 
di 128vrs 

/ cos 6 > b N 1 
V D v , 



cosOtfu sin6> t / 2 i 



Dr 
ti 



m 2 w -2(s + mZ ) 



-It \ S + t — ITlyy — 77j| J — 2m^m~ H — 



m 



w 



2cos^_ Re ( cosgt/n _ sme t f 21 ^ § , _ i + 2 ^ + 



w , 



t(js — m~ ^ 



Dr 



til 



W 



where Dr = s—mi +imr T~ and D v =t. In this formula, we have included diagrams 

not only for the t\ % exchanges in the s-channel but also for the neutrino exchange in the 
t-channel. 

In Fig. 2 we show the mr dependence of the total cross section for the process, 

e + p -> hW + X, (14) 



5 



where we take m~ = 205GeV, m~ = 225GeV and A' 131 = 0.1. Two sets of mixing angles 
(a) (8 t , 9b) = (0.95, 1.0) and (b) (1.2, 1.2) are adopted. They are suitable parameter 
sets for simulation of the excess events at the high Q 2 region in the NC DIS. Here we 
use MRS-G parton distribution [27]. We find that large cross sections of 0.5 ~ lpb are 
expected for 90GeV ^ m % 1 & 120GeV. The double threshold structure originated from 

both ti and ti contributions can be seen. 

There are two types of experimental signatures of the process (14) ; (A) a large missing 



transverse momentum fx plus multijets and (B) a high Pt charged lepton (e + or plus 
mono-jet with a large Ifx- (A) and (B) respectively correspond to hadronic decays and 
leptonic decays of the W + boson. In the next section, we will investigate the signal (A), 
which must be a clean signal of the process ( |i~I| ) because of Br (IF 4 " — > qq') = 2/3. It is 
worth mentioning that we expect a few events of the type (B) for an integrated luminosity 
of 30pb _1 and a typical efficiency of 50% due to a to t — lpb and Br(IF + — > vl + ) = 1/9. 



5 Event excess in CC process 

Now we consider one of signatures for the process, e + p — ► b\W + X, followed by the decays, 
b\ — > v e d and W + — > qq'. In this case a large Ifx plus multijets can emerge as its 
experimental signal. Such a sigal could be observed as an event excess at the high Q 2 
region in the CC DIS, because some multijets could not be distinguished from the mono- 
jet even after event selections of the CC DIS. In the event selection at HERA, a cut on 
Pt/Et is used to suppress backgrounds of the photoproduction and the NC DIS, where 
Pt = fx and Et = J2i E% sin#, (summed over the calorimeter cells). Indeed, lower cuts on 
Pt/Et is efficient in discriminating mono-jet events from multijet events b ecause Pt / Et 
~ 1 for the former case. However, the cuts on Pt/Et are not so severe in actual analyses at 
HERA, i.e., P T /E T > 0.4 and 0.5 at the ZEUS @ and HI |], respectively. Consequently, 
some mutijet events could be regarded as the CC DIS events after the selection. 

We investigate the Q 2 and M distributions of the expected events. In this analysis, 
we use the Jacquet-Blondel methods [K|, in which Q 2 B and Mjb are determined from 
energies and scattering angles of the final jets. We combine the HI and ZEUS data 
explicitly presented in refs.[16, 17|. Here we use MRS-G parton distribution [27] and the 
hadronization has been simulated a la JETSET [29]. 



Shown in Fig. 3 are Qjb and Mjb distributions of the expected number of events 
together with the experimental data. Figure 3 corresponds to the case that two stops 
are almost degenarate in mass (nij^, mj) = (205GeV, 225GeV) with the finite mixing 
angle 8t = 0.95. Parameters in the sbottom sector are taken as m~ = lOOGeV and Of, = 

1.0. In this case we obtain Br(ti -> e + d) ~ 20% and Br(t 2 -> e + d) ~ 10%. That is, the 
parameter set has not certainly been excluded by the leptoquark searches at Tevatron. 
Note, moreover, that the broad mass (or x) distribution of the excess events in the NC 
DIS 0, ^] can be reproduced by the mass spectrum and the mixing parameters ||. In 
Fig.3 we take Pt > 50GeV and Pt/Et > 0.5 as kinematical cuts. Additional cut, Q 2 > 
12,500GeV 2 , is adopted for Mjb distribution. The tendencies for event excess at the high 
Q 2 g and large Mjb of the data are successfully reproduced by our scenario though in the 
limited statistics. 
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We have to emphasize that the signal is characterized by the multijet events with 
a large fx- It means that the signal has a broad Pt/Et distribution. Obviously it is 
different from the one for CC DIS, which has a sharp peak at Pt/Et ^ 1. Therefore, 
we could extract the signal from the CC DIS by imposing a lower value of cut on Pt/Et 
such as Pt/Et > 0.3. Needless to say, detailed jet-analyses will also be indispensable for 
certainty. 

It should be mentioned that the possible event excesses in the Q 2 jb and Mjb distri- 
butions are originated from our choice of rather large coupling parameter A' 131 = 0.1. If 
more stringent upper limit, A' 131 & 0.05, will be established by elaborated analyses of the 



APV experiments 24], expected events from the signal would be slightly decreased. 



6 Concluding remarks 

We have investigated the process e + p — * b\W + X from the CC channel at HERA. Partic- 
ularly, we considered a case where the sbottom b is lighter than the stops and has a mass 
smaller than about 120GeV. Here we assume sufficiently heavy chargino W\ (<; 200GeV) 
and neutralino Z\ (<; lOOGeV) for simplicity. In this case the stops cannot decay into bW\ 
and the sbottom can only decay into i> e d via the R-parity breaking coupling A' 131 . We have 
shown that such a light sbottom has not been excluded by the present experiments. One 
of characteristic signatures of the process could be a large missing transverse momentum 
plus multijets corresponding to hadronic decays of the W . It could be observed as an 
event excess at the high Q 2 region in the CC DIS. We have shown that expected event 
distributions are consistent with the CC data recently obtained by the HI and ZEUS 
groups at HERA. Another signature of our process should be a high Pt charged lepton 
(e + or fi + ) plus mono-jet with a large ^t corresponding to leptonic decays of the W. We 
can expect a few events of such a characteristic signal assuming the integrated luminosity 
of 30pb _1 . Detailed studies of the leptonic signal of our process are in progress. 
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Figure 2: dependence of the total cross section a(e + p — ► biW+X). We take A' 131 

0.1, m T = 205GeV and m 7 = 225GeV. 
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Figure 3: Q JB and Mjb distribution of the expected number of events together with the 
experimental data. We take m z = 205GeV, m~ = 225GeV, = lOOGeV, 9 t = 0.95, 



In 



6b = 1.0, A' 131 = 0.1 and integrated luminosity L = 34pb . As kinematical cuts, Pt > 
50GeV and P T /E T > 0.5 are adopted. Additional cut, Q 2 > 12,500GeV 2 , is taken for 
Mjb distribution. Dashed line corresponds to the SM CC expectation. 
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